home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / program / nrpas13.zip / FLMOON.PAS < prev    next >
Pascal/Delphi Source File  |  1991-04-29  |  814b  |  25 lines
  1. PROCEDURE flmoon(n,nph: integer; VAR jd: integer; VAR frac: real);
  2. VAR
  3.    i: integer;
  4.    rad,xtra,t2,t,c,as,am: real;
  5. BEGIN
  6.    rad := 3.14159265/180.0;
  7.    c := n+nph/4.0;
  8.    t := c/1236.85;
  9.    t2 := sqr(t);
  10.    as := 359.2242+29.105356*c;
  11.    am := 306.0253+385.816918*c+0.010730*t2;
  12.    jd := 2415020+28*n+7*nph;
  13.    xtra := 0.75933+1.53058868*c+(1.178e-4-1.55e-7*t)*t2;
  14.    IF ((nph = 0) OR (nph = 2)) THEN BEGIN
  15.       xtra := xtra+(0.1734-3.93e-4*t)*sin(rad*as)-0.4068*sin(rad*am)
  16.    END ELSE IF ((nph = 1) OR (nph = 3)) THEN BEGIN
  17.       xtra := xtra+(0.1721-4.0e-4*t)*sin(rad*as)-0.6280*sin(rad*am)
  18.    END ELSE BEGIN
  19.       writeln('pause in FLMOON - nph is unknown.'); readln
  20.    END;
  21.    IF (xtra >= 0.0) THEN i := trunc(xtra) ELSE i := trunc(xtra-1.0);
  22.    jd := jd+i;
  23.    frac := xtra-i
  24. END;
  25.